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1.  Abstract 


The  results  of  this  research  centered  on  the  experimental  studies  of  a  single 
superconducting  persistent  current  qubit,  the  implementation  of  type-II  algorithms  using 
these  qubits,  and  the  proposal  for  adiabatic  quantum  computing  using  these  qubits. 

The  major  experimental  results  on  single  superconducting  persistent  current  qubits  have 
been  the  observation  of  the  quantum  energy  level  crossings  in  niobium  qubits,  and  the 
microwave  measurements  of  intra-well  relaxation  times. 

We  have  developed  two  implementation  methods  for  solving  the  one-dimensional 
diffusion  equation  with  a  type-II  algorithm.  In  the  first  method,  the  state  of  each  qubit  is 
set  by  the  local  magnetic  field  bias.  Although  this  initialization  method  has  the  advantage 
of  simplieity,  the  subsequent  unitary  collision  operations  demand  precise  timing.  The 
second  method  uses  nearly  identical  qubits  that  can  be  addressed  locally  at  the  node. 
Microwave  pulses  are  then  used  for  the  initialization,  and  the  unitary  transformation  is 
simplified  to  just  timed  free  propagation. 

A  scalable  architecture  for  an  adiabatic  quantum  computer  has  been  proposed  for 
superconducting  persistent  current  qubits  in  whieh  an  adiabatically  varying  magnetie 
field  is  applied  to  all  the  qubits  simultaneously. 
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2.  Summary  of  Results 


The  previous  work  has  been  focused  on  three  areas:  (1)  the  use  of  superconducting  qubits 
to  implement  type-II  quantum  computing  [1-6];  (2)  in  concert  with  the  DURINT  program 
(the  joint  AFOSR-ARDA  MURI)  on  type-I  quantum  computing,  the  study  of  single  and 
coupled  superconducting  persistent  current  qubits[7,8];  and  (3)  the  proposal  of 
architectures  for  adiabatic  quantum  computing  [9]  using  superconducting  qubits. 

2.1  Persistent  Current  Qubit  (PC-  Qubit) 

The  particular  device  that  we  have  studied,  co-invented  by  one  of  us  [8,  10],  is  made 
from  a  loop  of  Nb  interrupted  by  3  Josephson  junctions  (Fig.  1).  The  application  of  an 
external  magnetic  field  to  the  loop  induces  a  circulating  current  whose  magnetic  field 
either  adds  to  (say  circulating  current  in  the  clockwise  direction)  or  opposes 
(counterclockwise)  the  applied  magnetic  field.  When  the  applied  field  is  near  to  one-half 
of  a  flux  quantum,  both  the  clockwise  and  counterclockwise  current  states  are  classically 
stable.  The  system  behaves  as  a  two-state  system.  The  samples  are  fabricated  at  MIT's 
Lincoln  Laboratory  in  niobium  by  photolithographic  techniques  on  a  trilayer  of  niobium- 
aluminum  oxide-niobium  [11]. 


I — X — ' 

' — X — I 

FIGURE  1.  (a)  SEM  image  of  the  persistent  current  qubit  (inner  loop)  surrounded  by  the  measuring 
dc  SQUID,  (b)  a  schematic  of  the  qubit  and  measuring  SQUID,  the  x's  mark  the  Josephson 
junctions,  (c)  The  energy  levels  for  the  ground  state  (dark  line)  and  the  first  excited  state  of  the  qubit 
versus  applied  flux.  The  double  well  potentials  are  shown  schematically  above.  The  lower  graph 
shows  the  circulating  current  in  the  qubit  for  both  states  as  a  function  of  applied  flux.  The  units  of 
flux  are  given  in  terms  of  the  flux  quantum. 
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The  Hamiltonian  for  the  single  PC  qubit  is 


^pcubit 

where  /p  is  the  circulating  current  in  the  qubit  and  /  is  the  magnetic  flux  threading  the 
loop. 


Our  work  on  single  qubits  leverages  off  our  work  on  type-I  quantum  computing,  and  our 

major  experimental  results  have  been 

a.  Experimental  mapping  of  the  potential  energy  landscape  of  the  persistent  current 
qubit  by  thermal  activation  studies[  12]. 

b.  Observation  of  the  quantum  energy  level  crossings  at  20  mK  in  niobium  qubits  fabri¬ 
cated  at  MIT  Lincoln  Laboratory. 

c.  Microwave  spectroscopy  of  the  energy  levels  in  the  qubits,  which  agrees  with  the 
quantum  energy  level  calculations.  Long  intra-well  relaxation  times  greater  than  50 
microseconds  were  observed,  as  well  as  and  T2*  times  of  about  10  nanoseconds.  Both 
of  these  times  will  improve  with  improved  environmental  protection. 

2.2  Type-II  quantum  computing 

We  have  developed  two  implementation  methods  for  solving  the  one-dimensional 

diffusion  equation  with  a  type-II  algorithm: 

1.  Each  type-II  algorithm  vmdergoes  the  following  stages:  (1)  the  initialization  of  the 
coupled  qubits  in  a  node;  (2)  the  collision  operator  which  performs  a  quantum 
unitary  transformation  at  the  node;  (3)  the  measurement  of  the  qubits;  and  (4)  the 
streaming  of  classical  information.  The  hallmark  of  die  first  method  that  we 
developed  is  that  the  state  of  each  qubit  is  initialized  by  a  static,  local  magnetic 
field  bias,.  Although  this  initialization  method  has  the  advantage  of  simplicity,  the 
subsequent  unitary  collision  operations  demand  precise  timing. 

2.  The  second  method  maps  our  superconducting  system  onto  tire  NMR  system  tiiat 
was  used  to  demonstrate  the  type-II  algorithm  for  the  diffusion  equation  [13]. 
This  mediod  uses  nearly  identical  qubits  diat  can  be  addressed  locally  at  the  node. 
Microwave  pulses  are  then  used  for  the  initialization,  and  the  unitary 
transformation  is  simplified  to  just  timed  free  propagation.  Circuits  of  two 
coupled  qubits  have  been  designed  and  are  being  febricated  to  study  this  method. 


a.  The  algorithm 

The  simplest  type  n  algorithm  to  implement  is  the  Factorized  Quantum  Lattice-Gas 
Algorithm  (FQLGA)  for  the  ID  Diffusion  Equation  [1,4].  It  requires  a  ID  network  of 
classically  coimected  quantum  computers  (nodes)  that  consist  of  only  two  coupled  qubits. 
Now  each  node  represents  a  position  in  a  ID  space,  with  the  four  quantum  states  at  each 
spatial  position  representing  a  completely  unoccupied  spatial  position  (|0)a|0)b),  those 
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with  a  particle  moving  to  the  left  only  (|l)a|0)b),  those  with  a  particle  moving  to  the  right 
only  (|0)a|l>b),  and  those  with  particles  moving  in  both  directions  (|l)a|l)b).  The  only 
collision  that  can  occur  between  “particles”  is  for  the  scenario  where  the  amplitudes  of 
the  |0>a|l)b  and  the  |l>a|0)b  states  are  exchanged.  The  collision  thus  simply  swaps  the 
middle  two  particle  occupation  scenarios.  The  collision  operator  that  implements  this 
exchange  of  amplitudes  we  will  call  sqrtswap,  which  will  be  described  in  more  detail 
below.  As  for  the  number  density  of  the  fluid,  it  is  just  the  sum  of  the  probability  for 
each  qubit  to  be  in  the  occupied  state. 

Note  that  it  is  not  obvious  that  the  prescription  just  described  should  reproduce  ID 
Diffusion  equation.  There  are  lattice  requirements  and  spurious  invariant  requirements 
not  discussed  here  which  must  be  satisfied  as  well. 

We  will  now  go  through  each  step  in  the  algorithm  in  detail. 

The  first  step  of  the  algorithm  is  the  initialization  portion,  which  sets  the  pre-collision 
number  density  p  =  Pa  +  Pb,  where  P  is  the  probability  for  a  qubit  to  be  in  the  excited  i.e. 
occupied  state: 


r,=>/r^,|o)+v^|i) 

Qubit  A  represents  the  occupancy  of  a  particle  moving  to  the  left  and  qubit  B  represents 
the  occupancy  of  a  particle  moving  to  the  right. 
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FIGURE  2.  Schematic  of  sequential  operations  for  type-II  quantum  computation  algorithms 

The  algorithm  steps  can  be  visualized  as  shown  above.  The  first  step  is  to  initialize  each 
of  the  qubits  in  each  node.  For  simplicity,  there  are  6  nodes  shown,  and  each  node  has 
two  qubits.  The  collision  portion  of  the  fluid  dynamics  is  represented  by  the  sqrtswap 
imitary  transformation  written  as: 


2  0  0  0^ 

0  1+/  \-i  0 
0  1  —  z  1  +  z  0 
0  0  0  2; 

After  the  collision  the  system  can  be  represented  as  being  in  a  single  state.  The  third  step 
is  to  measure  the  state  of  each  qubit  and  then  repeat  the  above  three  steps  to  produce  a 
mean  occupancy  number  P’  for  every  qubit.  After  the  calculation  of  P’  is  complete,  the 
density  of  the  fluid  can  be  determined  by  summing  the  mean  occupancies  of  all  Ae  qubits 
at  each  node.  Lastly,  repeat  the  above  for  the  next  time  step,  initializing  Ihe  number 
density  with  that  from  the  appropriate  moving  particle  from  the  previous  time  step.  Two 
methods  to  implement  this  series  of  steps  for  the  ID  Diffusion  equation  are  discussed 
below. 
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b.  First  Implementation  Concept 

Without  considering  the  common  initialization  and  transformation  concepts  of  quantum 
computing,  this  algorithm  seems  like  it  could  be  implemented  in  a  rather  simple  way. 
Ignoring  all  classical  necessities  like  storing  the  individual  measurement  results, 
computing  the  average  of  many  of  these,  streaming  the  occupation  niimbers,  etc.,  the  only 
quantum  mechanical  steps  are  the  initialization,  the  collision,  and  the  measurement. 

Let  us  begin  with  the  simplest  step,  the  measurement.  This  can  be  accomplished  with 
an  underdamped  DC  SQUID  with  our  typical  ramping/switching  measurement  technique. 
The  SQUID  can  distinguish  the  two  states  of  opposite  circulating  current  due  to  die 
opposite  fluxes  they  produce,  hence  distinguishing  CTz  eigenstates. 

The  initialization  portion  of  the  algorithm  can  be  done  by  taking  advantage  of  the 
constant  Ox  term  in  our  Hamiltonian,  which  is  added  to  the  familiar  NMR  cTz  term  that 
depends  linearly  on  your  tunable  external  field,  here  being  the  flux  (J)  within  the  qubit  as 
opposed  to  the  applied  magnetic  field  in  NMR.  Initialization  is  then  accomplished  by 
letting  die  qubit  relax  to  its  ground  state  and  setting  the  flux  which  in  turn  sets  the  qubit 
in  any  of  the  continuum  of  eigenstates  that  range  from  |0)  to  |1)  along  the  Bloch  sphere 
geodesic  with  real  phase. 
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FIGURE  3.  (a)  Typical  avoided  crossing  diagram  showing  two  bands  and  the  ground  and  excited 
states  in  terms  of  the  asymptote  states;  (b)  Occupancy  of  ground  state  as  a  function  of/. 

This  type  of  initialization  is  done  while  the  qubits  are  uncoupled,  which  can  either  be 
remedied  by  actually  decoupling  the  qubits  during  initialization,  or  by  just  initializing 
them  while  coupled,  which  would  be  feasible.  ^ 

The  collision  at  first  glance  looks  very  simple  too,  but  on  closer  consideration  proves  w 
be  demanding.  The  sqrtswap  transformation  simply  swaps  the  states  |0)a|l)b  and  \\)^^, 
i.e.  performs  a  Rabi  oscillation  between  these  two  states  such  that  ORabit  =  rt.  The  first 
problem  arises  because  our  initialization  technique  has  set  our  coupled  Hamiltonian  such 
that  the  middle  two  eigenstates  in  energy  are  most  likely  not  |0)a|l)b  and  |l)a|0)b,  therefore 


8 


applying  radiation  at  the  frequency  corresponding  to  the  middle  eigenstates’  energy 
splitting  for  ©Rabi^  =  will  not  perform  the  desired  transformation.  One  way  around  this 
is  to  find  an  operating  point  it,)  where  the  middle  two  eigenstates  are  veiy  nearly 
|0>a|l)b  and  |l)a|0>b,  such  that  applying  radiation  at  their  energy  difference  for  t  =  7r/©rabi 
will  accomplish  the  collision.  Below  is  one  example  of  an  operating  point  where  the  two 
sets  of  eigenstates  are  very  similar;  there  |2)  is  the  first  excited  state  of  the  coupled 
system  and  |3)  is  the  second  excited  state  of  the  coupled  system. 


FIGURE  4.  Overlap  of  middle  two  eigenstates  and  middle  two  computational  basis  states 

We  have  verified  both  that  the  energy  difference  between  these  two  states  at  this 
operating  point  as  well  as  the  matrix  elements  connecting  the  middle  two  eigenstates  are 
non-zero,  so  ever3dhing  is  fine  so  far.  We  have  also  verified  that  fixe  error  introduced 
from  this  approximation  of  the  collision  does  not  undermine  the  results  of  fluid 
simulation  altogether  as  can  be  seen  below.  We  have  plotted  the  number  density  p  on  the 
ordinate. 


FIGURE  5.  Comparison  of  swap  operations  and  influence  on  type-II  QC 

Once  the  initial  states  are  set  (with  the  qubits  uncoupled),  the  subsequent  radiation  has 
to  be  turned  on  precisely  when  you  couple  the  qubits  since  the  qubit  states  will  begin 
evolving  immediately  upon  coupling,  on  larmor  time  scales.  Dealing  with  larmor  time- 
scale  timing  is  experimentally  demanding,  making  this  option  very  unappealing. 

One  could  also  initialize  the  qubits  while  they  are  coupled,  with  only  a  small  amount  of 
error,  and  remove  the  above  constraint,  but  there  is  another  constraint  that  appears  in 
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either  initialization  technique.  To  change  the  operating  points  from  the  initializing  values 
to  the  values  that  will  produce  a  good  collision,  one  must  start  the  radiation  immediately 
after  this  is  done,  on  a  larmor  time  scale.  To  avoid  having  to  control  our  qubit  on  this 
rapid  time  scale,  we  have  developed  an  initialization  and  control  method  similar  to  that 
used  in  NMR. 


c.  NMR-like  Implementation  Concept 

Initialization  in  NMR  is  typically  done  by  applying  radiation  at  the  frequency 
corresponding  to  the  energy  difference  between  the  two  single  qubit  states,  with  the 
coupling  term  ignored  as  long  as  the  radiation  amplitude  (expressed  as  an  energy)  is 
sizable  compared  to  the  coupling  constant.  In  the  rotating  frame  that  is  spinning  at  the 
radiation  frequency  around  the  z-axis,  with  radiation  applitd  on  resonance,  the  NMR 
Hamiltonian  is: 


H  nmr  —  iYfi[cos(</5p)4  +  sin(^p)/y] 

JU 

where  ^  is  the  phase  of  the  radiation,  y  is  the  magnitude  of  the  nuclei’s  magnetic 
moment,  and  B  is  the  magnitude  of  the  applied  radio  frequency  magnetic  field. 

With  this  Hamiltonian  you  can  now  perform  arbitrary  single  qubit  rotations  in  the 
rotating  frame  by  applying  radiation  for  times  that  yield  your  angle  of  rotation  (0)  via 


using  different  initial  phases  to  choose  the  axis  of  rotation. 

To  do  a  sequence  of  two  orthogonal  axis  rotations  for  different  times  to  initialize  the 
qubit,  one  need  not  have  larmor  time  scale  precision.  For  the  three  junction  persistent 
current  qubit  applying  radiation  on  resonance,  in  the  rotating  frame  spinning  at  the 
radiation  frequency  around  the  stationary  eigenvalue  axis  z’  gives, 

/fpcqubit  =  —  gQsin(^)[cos(^p)/,(.  +sin(^p)/y] 

where  go  is  proportional  to  the  amplitude  of  your  radiation  and  0  defines  your  operating 
point.  Atf=<I>o/2, 

=  11 

H pcqubitCy 

where  now  your  rotating  frame  is  spinning  around  the  x-axis. 

Measurement  of  the  qubit  state  must  be  performed  in  the  Ox  basis  so  that  all 
computation  can  be  done  in  the  rotating  frame  atf=  0(/2.  Operating  at^Oo/2  is  ideal 
for  low  flux  noise,  but  some  care  must  be  taken  about  the  measurement.  At  first  glance  it 
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would  seem  that  the  eigenstates  are  not  resolvable  with  the  measurement,  but  due  to 
self-inductance  effects  this  is  not  true. 

With  the  radiation  off,  in  the  lab  frame,  the  Hamiltonian  for  a  two  qubit  NMR  system 
reads: 


Hmm  = 

Since  the  measurement  is  in  the  Oz  basis,  and  to  preserve  computing  in  the  rotating  frame, 
you  always  go  to  a  rotating  frame  about  the  z-axes.  In  a  co-rotating  frame  at  the  resonant 
frequency: 


n 


From  here  it  is  straightforward  to  get  to  any  other  imaginable  two  qubit  operation. 
For  the  persistent  current  qubit  things  are  a  bit  more  complicated  since  the  coupling  term 
is  not  in  the  same  direction  as  the  single  qubit  terms  at/=  <l>o/2  : 


H  pcqubit 


iTt 


J  PP 


Preserving  the  ability  to  compute  in  the  rotating  frame  forces  one  to  rotate  both  qubits 
about  the  x-axis  at/=  3)o/2,  making  the  transformation  less  trivial  than  in  the  NMR  case, 
resulting  in: 


H  pcqubit  —~J\2  U z^ z  1 

n 


where  the  constraint  that  coi  =  (02  has  been  imposed.  This  removes  the  ability  to  address 
the  qubits  separately  via  frequency  discrimination,  and  forces  one  to  address  them 
separately  spatially,  i.e.  each  seeing  a  different  radiation  source.  Fortunately,  one  benefit 
of  solid  state  computing  is  that  each  qubit  can  have  its  own  personal  antenna. 

Now  as  to  the  exact  collision  operator  we  need  to  implement,  we  have  the  following: 


1 

sqrtswap  =  — 


{2  0  0  0^ 
0  l+i  l-i  0 
0  l-i  1  -t-  i  0 
0  0  0  2 


=  exp[-i-|(cr’o-,^  -f-  cr^cr^  -1-  a^arl )] 

O 


Since  the  arguments  in  the  exponential  commute  you  can  write  the  collision  matrix  as: 


sqrtswap  =  -i-  o-^o-^)]exp[-i-J<T>^^] 

O  O 

The  first  term  is  just  free  evolution  in  the  doubly  rotating  frame.  The  second  term  can  be 
written  as: 
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The  middle  term  on  the  right  hand  side  can  then  be  achieved  in  the  rotating  frame  by 


expH^cr^crf  ]  =  expH-^(<T^crf  +cr’cry)]/?2(^)expH^(<j‘crf  +  cr\<7l)'\R],{n:) 
8  O  O 


And  so  at  the  price  of  locally  addressing  each  qubit  with  radiation,  which  was  done  in  the 
first  scheme  except  with  a  static  field  for  initialization,  now  only  pulses  on  the  Rabi  time 
scale  are  needed,  and  none  on  the  larmor  scale. 


In  summary,  the  following  design  criteria  must  be  met  for  the  second  implementation 
scheme: 


1.  The  two  qubits  must  be  identical  so  that  your  coupled  free  evolution  is  static  in 
the  doubly  rotating  frame. 

2.  The  radiation  must  be  applied  such  that  it  is  mostly  seen  by  only  the  desired  qubit, 
one  radiation  line  for  each  qubit. 

3.  The  radiation  must  be  at  least  50  times  greater  in  magnitude  than  the  coupling 
constant,  so  that  the  coupled  static  term  in  die  rotating  frame  can  be  ignored 
during  initialization  and  any  other  single  qubit  operations. 

4.  The  DC  SQUID  must  have  no  resonances  at  the  energy  splitting  of  the  qubits. 


2.3.  Adiabatic  quantum  computing 

A  scalable  architecture  for  an  adiabatic  quantum  computer  has  been  proposed  for 
superconducting  persistent  current  qubits  in  which  the  computation  is  performed  by 
slowly  varying  a  magnetic  field  applied  to  all  the  qubits  simultaneously[14].  Only  local 
incoherent  operations  are  needed,  such  as  switching  on  and  off  the  pairwise  interaction, 
and  measurement  of  only  a  small  subset  of  the  qubits. 

Adiabatic  quantum  computation  [9]  is  a  recently  proposed,  general  approach  to  solving 
computational  problems  of  the  complexity  class  NP  via  energy  minimization.  In 
particular,  by  exploiting  the  ability  of  coherent  quantum  systems  to  follow  adiabatically 
the  ground  state  of  a  slowly  changing  Hamiltonian,  it  aims  to  bypass  automatically  the 
many  separated  local  minima  that  occur  in  difficult  minimization  problems  and  confound 
all  known  classical  heuristics.  Adiabatic  quantum  computation  is  of  theoretical  interest 
because  it  provides  a  straightforward,  non-oracular  way  to  pose  class  NP  problems  on  a 
quantum  computer,  and  most  research  on  it  to  date  focuses  on  ascertaining  its  time 
complexity  [9,  15,  16].  However,  it  is  also  of  potentially  great  practical  interest  because 
encoding  a  quantum  computation  in  a  single  eigenstate,  the  ground  state,  offers  intrinsic 
protection  against  dephasing  and  dissipation  [17]. 

We  also  discuss  the  implementation  scheme  of  Kaminsky,  Lloyd,  and  Orlando  for 
using  superconducting  persistent-flux  qubit.  This  proposed  architecture  is  robust  in  the 
face  of  manufacturing  imperfections,  decoherence,  and  noise.  They  found  the  following 
criteria  for  adiabatic  QC  architecture: 

1 .  A  starting  Hamiltonian  with  an  easily  initialized  ground  state. 
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2.  Initial  and  problem  Hamiltonians  that  share  no  symmetries  so  as  to  ensure  levels 
do  not  cross. 

3.  To  minimize  unwanted  components  to  the  interqubit  couplings 

4.  To  equalize  the  mutual  inductances  between  qubits  to  be  coupled,  even  if  they  are 
not  nearest  neighbors 

5.  To  allow  room  for  the  shrinkage  of  the  energy  gap  expected  when  interpolating 
between  Hamiltonians. 

To  have  a  starting  Hamiltonian  with  an  easily  initialized  ground  state,  initialize  the 
qubits  with  f  =  1/2  flux  quanta  so  that  each  qubit’s  term  dominates  all  other  terms 
affecting  the  qubit.  This  choice  allows  the  interqubit  couplings  to  be  switched  on  for  the 
entire  computation.  Switching  couplings  on  and  off  is  needed  only  to  program  the 
computer;  it  is  not  neeessary  during  any  computation. 

To  have  the  initial  and  problem  Hamiltonians  that  share  no  symmetries  so  as  to  ensure 
levels  do  not  cross,  use  a  4-junction,  rather  than  3-junction  PC-qubit.  The  four-junction 
qubit  replaces  the  qubit  with  the  lower  critical  current  with  a  tunable  parallel  junction,  so 
as  to  allow  XX  and  XZ  inductive  couplings  as  well  as  ZZ. 

To  minimize  unwanted  eomponents  to  the  interqubit  couplings,  use  +  <7x  as  basis  of 
the  problem  Hamiltonian.  This  naturally  incorporates  the  inevitable  mixed  02  ox  terms 
when  4-junction  qubits  are  inductively  coupled  simply  by  having  the  coupling  loop 
arranged  symmetrically  so  as  to  have  equal  mutual  inductances  to  both  loops  of  the  4- 
junction  qubit,  since; 

(ozi  +  crxi)(crz2  +  0x2)  =  ctziOzi  +  0210x2  +  OxiOz2  +  0x10x2 


FIGURE  5.  schematic:  not  to  scale,  inductive  loop  much  longer  than  qubit  dimensions 


To  equalize  the  mutual  inductances  between  qubits  to  be  coupled,  even  if  they  are  not 
nearest  neighbors,  it  is  essential  that  the  buses  through  which  non-nearest  neighbors  are 
coupled  have  as  low  self-inductances  as  possible.  This  criteria  follows  from  the  fact  that 
the  qubit-qubit  mutual  inductance  via  a  coupling  loop  goes  ashe  ratio  of  the  square  of  the 
mutual  inductance  of  the  qubit-coupling  loop  to  the  self  inductance  of  the  coupling 
loop. A  combination  of  minimizing  the  loop  area  of  the  buses,  twisting  their  wires,  and/or 
placing  ground  planes  beneath  them  should  accomplish  this. 


Lastly,  to  allow  room  for  the  shrinkage  of  the  energy  gap  expected  when  interpolating 
between  Hamiltonians,  Optimize  junction  parameters  so  as  to  maximize  tunnel  splitting 
and  critical  current  (so  as  to  maximize  Mf)  .  For  example,  numerical  calculations  on  the 
full  4-junction  qubit  Hamiltonian  for  our  current  technology  of  Nb  junctions  of  Jc  =  2 
(iA/pm^of  sizes  (0.85  pm)2  and  (0.65  pm)2  yields: 

( 
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Energy  Difference  between  Qubit  Levels  =  3.7  GHz  (190  mK) 

Inductive  Coupling  »  4.5  GHz/pH  (225  mK/pH) 

Thus,  given  the  apparent  behavior  of  gap,  current  technology  should  permit 

demonstrations  of  the  ideal  adiabatic  algorithm  with  several  qubits  and  of  the  tunneling- 
based  heuristic  with  many  more. 

Future  Josephson  junction  technology  could  provide  Josephson  energies  10-100  times 
higher  than  those  currently  used  in  our  PC  qubits  and  thus  allow  for  ideal  adiabatic 
computers  with  lO’s  or  lOO’s  of  qubits,  and  quantum  annealers  with  lOOO’s  of  qubits. 
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